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Using Finite-Difference Time-Domain (FDTD) simulation, we show that ultrahigh- 
Q nanocavities can be obtained through the manipulation of a single semiconductor 

O. 

nanowire (NW) inside a slot in a line defect of a two-dimensional (2D) photonic 



crystal. By controlling the design and its lattice parameters of the photonic crystal, 



we have achieved a quality factor Q larger than 10 6 and a mode volume V c smaller 



than 0.11 /iin 3 (1.25 of a cubic wavelength in the NW) for a cavity peak in the 
telecommunication band. This design is useful for realizing a position-controlled 
cavity in a photonic crystal. Here we also discuss the small dependence of the Q- 
factor, the V c , and the cavity peak in relation to the position of the NW inside the slot 
and the potential application to the cavity quantum electrodynamics (QED) using 
the embedded-emitter NW. 
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I. INTRODUCTION 



In recent years, many new types of optical cavities with a high quality factor (Q) and a 
small mode volume (V c ) have been developed and fabricated^. From among these designs, 
photonic crystal slab nanocavities have attracted a lot of attention since they offer a very 
high Q/V c and the capacity for integration with other devices through nanowaveguides 2 . 
These nanocavities may also exhibit an ultrahigh-Q(> 10 6 ) and a small (A/n) 3 )^-, and 
this has led to the proposal and demonstration of ultralow power devices including power 
switches^, ultralow threshold lasers^, and all-optical random access memory (RAM) 1 ^. 
Thus, the high-Q nanocavities may enable the dense integration of ultralow power photonic 
devices in a single chipii. 

Most ultrahigh-Q nanocavities are fabricated by a combination of high-resolution electron 
beam lithography and dry etching and the nanocavities must be incorporated in the photonic 
crystal design prior the fabrication process. But recently, the high-Q nanocavities have been 
obtained through a post modification of photonic crystal waveguides by employing scanning 
probe lithography 12 . Some high-Q (< 10 4 ) nanocavities can also be created by filling the 
holes in the photonic crystals via water or polymer micro-infiltration 1 ^ 4 - but the feasibility 
of realizing a liquid-filled cavity is relatively low. Also, although these methods provide a 
rewritable cavity, there is no method with the ability to control the position of the cavity. 
Such a movable cavity would be useful in on-demand position-controlled devices in photonic 
integrated circuits. 

A single semiconductor nanowire (NW) in the line defect of a two-dimensional (2D) 
photonic crystal may offer the possibility of realizing a position-controlled cavity since a 
single NW may modify the refractive index at any position in the line defect. The coupling of 
a single semiconductor NW embedded in a one-dimensional (ID) photonic crystal is intended 
solely to reduce the cavity losses of the NW caused by the reduction of the evanescent 
field outside the NW body^^. Furthermore, the coupling of an NW with a 2D photonic 
crystal waveguide is designed solely for efficient photon delivery 1 ^- 1 ^. For coupling between 
a cavity and an emitter, NWs with an embedded single emitter (e. g. quantum dots 
(QDs) 1 ^— and diamond color centers 2 ^) already exhibit excellent performances as single 
photon sources. Therefore, in addition to the position-controlled cavity, the incorporation 
of a single embedded emitter in a single NW may also provide a site-controlled emitter for 
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a light matter interaction. This position-controlled cavity and emitter could be interesting 
as regards the revolution taking place in the fields of quantum computation^, quantum 
cryptography 27 , and quantum optics^. 

In this work, we propose high-Q nanocavities based on a single semiconductor NW inside 
the air slot of a line defect in a 2D photonic crystal, and we demonstrate nanocavities with 
Q ~ 2 x 10 6 and V c < 0.11/zm 3 (1.25 of a cubic wavelength in the NW) for the cavity 
peak in the telecommunication band. Our approach shows the possibility of engineering 
the position of the NW through NW manipulations. We use a numerical modeling with 
the three-dimensional Finite-Difference Time-Domain (3D FDTD) method to show that our 
platform still exhibits a relatively large Q and small V c despite abrupt changes in the NW 
size, slot depth, slot width, and NW position. We also estimate that the embedded QD 
NWs are well suited for cavity quantum electrodynamics (QED) experiments by considering 
the light matter interaction with the single QD at the center of the NW. 

II. CAVITY DESIGN 

In this work, we consider NWs with a refractive index (n^w) of 3.17 (i.e. InP NWs), 
while the substrate has a refractive index (n su b s ) of 3.48 (i.e. silicon (Si)). However, this 
approach is suitable for different material systems using different parameters. We assume a 
square cross-section for our NWs with a side length (si) of 50 nm and lengths (L) between 
3 and 5 /xm (see Fig. [T^i). These are conventional parameters for InP NWs grown using 
low-pressure metalorganic vapor phase epitaxy (MOVPE) with [001] direction on InP (001) 
substrate^. We choose a NW with a square cross-section to simplify the geometry analysis 
needed for our calculations. In Fig. [TJi shows a photonic crystal design based on a line-defect 
waveguide (consisting of a row of missing holes in the Y — K direction). The hole radius 
(R) and the thickness (t) are constants of 130 and 150 nm, respectively, while the lattice 
constant (a) is 380 nm. The width of the line defect (the distance between the nearest 
neighbor holes on either side of the missing-hole row) is 0.98xa\/3- For the air slot, we vary 
the parameters such as the slot depth (d = 25-75 nm), width (w = 50-200 nm), and the 
position of the NW to investigate the change in Q and V c caused by errors in slot fabrication 
or the nanopositioning of the NW. 

In our design, mode-gap confined photonic crystal nanocavities are realized by the NW 
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while the air slot is introduced as a dielectric discontinuity to enhance the electric field 
from the original cavities. The physical mechanism is very similar to that reported in 
previous studies on entirely-drilled air slots22^, in which the local field enhancement at 
dielectric boundaries is considered. However, unlike the cavity with one-sided cladding^, 
the termination of the air slot does not cause much Q degradation in spite of the strong 
vertical asymmetry in a 2D photonic crystal slab. The radiation loss to the lower part of 
the slab is suppressed by the NW in the center of the air slot, which increases the effective 
refractive index. 

In general, the electric and magnetic fields of the transverse electric (TE) mode in a 
photonic crystal cavity are mainly x— and y—, z— polarized, respectively. First, we consider 
electric field enhancement for the simplest photonic crystal design, where the NW side length 
matches the air-slot width. For maximum enhancement, w should be infinitesimally narrow 
while d should be shallow enough. It also means that the NW should be as thin as possible. 
A large slot and NW may reduce the field strength at the slot boundaries. This explains the 
choice of a NW with a 50-nm long side NW in our design since this size is suitably small 
for the NW. In the real design, the NW can not fit with the air slot perfectly. Then, the 
enhancement condition will be different in this case. For this reason, the w and d values of 
the air slot in our design must be optimized with the si of the NW so that the electric field 
of the cavity can be enhanced optimally. We also note that the misalignment of the NW 
from the center of the slot may reduce the field strength. 

The presence of the air slot in our design has certain advantages regarding the difficulty 
involved in the high precision nanopositioning of the NW in the line defect. Placing a 
single nanowire in the center of a simple line defect in the photonic crystal can create a 
cavity but with a moderate Q < 10 4 . However, the Q is extremely sensitive with respect 
to the position of the NW, especially the x-direction(see Fig. [CJi). The exact position 
of the NW in the center of the line defect requires more precision than that provided by 
current nanopositioning technology. Therefore, the air slot acts as a barrier for the spatial 
movement of the NW. With the air slot, we can also move the NW in the line defect along 
the z-direction with a small error in the x-direction. 

We can also control the refractive index of the line defect waveguide through the depth of 
the air slot. Theoretically, ultrahigh-Q nanocavity modes can be created by the ultrasmall 
local modification of the refractive index of a line defect waveguide in a 2D photonic crystal 
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slab 33 . This was demonstrated by employing a photosensitive process for chalcogenide glass 
photonic crystals 34 and scanning probe lithography for Si photonic crystals^. Here we 
propose that removing small amounts of Si in the line defect through the air slot affects the 
confinement around the NW. This mechanism is explained by the schematic band diagram in 
Fig. [Tb, which was calculated by the plane-wave expansion method. For this band diagram 
calculation, we set a = 380 nm, L = 3 /im, R = 130 nm, and si = w = d = 50 nm. Without 
the NW and the air slot, the line defect in the photonic crystal has a sharp mode gap edge 
in the band diagram (see black dotted lines in Fig. [Tb). If we lower the mode-gap edge of 
the central part through the air slot and the NW (red dotted lines), we create a confined 
mode surrounded by the mode gap in an unmodified line defect (green dashed lines). Such 
modification is equivalent to structural 3 -^ and refractive-index modulations 3 - 3 - 1 ^. 

To investigate the field distributions, we performed a 3D FDTD with a calculation area 
as large as 32.4 (x axis), 7.9 (y axis) and 40.3 (z axis) lattice periods. The perfectly matched 
layer method with 8 layers was employed as absorbed boundary condition in all directions. 
We set a grid spacing of 25 nm so that one of more field components was contained within 
the 50-nm NW side length. For the details of the method, we first performed the simulation 
using a broadband excitation light source and obtained a cavity spectrum through Fourier- 
transformed field. We defined the resonance angular frequency u and used this parameter 
to simulate the field and the energy decay of the cavity mode with a narrow band excitation. 
To determine Q which represents the strength of the light confinement, we fitted the tail 
of the energy decay curve with a single exponential and we determined the cavity photon 
lifetime r p h- We estimated the Q in a cavity as Q = r p h ■ u. 

III. FDTD SIMULATION FOR NW IN WIDTH-MATCHED AIR SLOT 

First, we simulate the field distributions for a condition where the side length of the NW 
perfectly matches the air-slot width in the line defect of the photonic crystal si = w = d = 
50 nm. Fig. [2K shows the electric- field {\E\) and magnetic- field {\H\) distributions in the 
x-z plane for the fundamental TE mode. The electric- and the magnetic-fields are mainly 
dominated by x- and y-polarizations, respectively, and both fields are confined in the 3-/zm- 
length NW. For this cavity mode, the cavity wavelength is exist at the telecommunication 
wavelength (1410 nm) while Q and V c are 4 x 10 5 and 0.07 /im 3 , respectively. Although 
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this design already exhibits high Q, an NW-width-matched air slot is difficult to realize 
experimentally, i.e. nanomanipulation of the NW requires a slot, that is larger than the 
NW. 



IV. FDTD SIMULATION FOR THE NW WITH DIFFERENT LENGTHS L 
IN THE WIDTH-MATCHED AIR SLOT 

It is important to simulate the sensitivity of the L value of the NW with Q and V c . 
Principally, the NW length can be determined by controlling the growing condition but it is 
difficult to choose a precise length of 3-/zm from many NWs with certain length distributions. 
Here we discuss the fidelity of Q and V c with different L values for our design from an analysis 
of the z-component of the magnetic-field (H z ) distributions in Fig. [3k- When the NW 
becomes longer, the mode distributions become spread along the waveguide (z) direction. 
We found that the full widths at half maximum (FWHM) of the mode distributions were 
1.46, 1.81, and 2.33 /im for 3-, 4-, and 5-/im NW, respectively. This shows that the 
spread of the mode distributions is proportional to L. From the profile along the waveguide 
depicted in Fig. [3b, as L increases, the intensity in the center decreases while that at the side 
increases. Quantitatively, Q/V c is relatively constant, which means that the enhancement is 
not sensitive to L (see Fig. |3b). The Q itself has a small variation. This is most likely due to 
the small difference in the local index modulation caused by different L in comparison with 
the large difference in the refractive index of NW and that of the photonic crystal. We 
observed that the increase in V c is accompanied by the spreading of the mode distributions, 
which is similar with the case of the local index modulation in the waveguide of the photonic 
crystal^. 



V. FDTD SIMULATION FOR NW WITH DIFFERENT SLOT DEPTHS d 

It is more difficult to fabricate an air slot with a partial slab depth (d) than an entirely- 
drilled air slot. An electron beam lithography or dry etching process may provide a different 
d from that expected. Therefore, we include the d variation (rf±0.5sZ) in our simulation for 
the investigation of Q and V c . In this calculation, other than d, we use the same parameters 
as in Fig. [2k,. 

6 



The magnetic-field distributions (H y ) in the x — y plane are shown in H^t. From these 
results, we observe that the field profiles along the ^/-direction are relatively similar while 
the profiles along the x-direction change with the position of the NW (see Fig. Hb and c 
for y- and ^-directions, respectively). The mode distributions have the longest spread in 
the x-direction with the shallowest air slot. The profile becomes shorter for a deeper air 
slot. With both shallow and deep air slots, the intensity in the center decreases. The cavity 
peak, Q and V c are summarized in Fig. HH. We notice that the cavity peak A exhibited an 
8-nm-red shift with a change in d. Therefore, this design is interesting in terms of detuning 
the wavelength using an air slot whose depth varies along the ^-direction. V c is relatively 
constant for different d values. The highest Q is found when the depth is the same as the 
NW side length d = si. The maximum decrease in Q is only by a factor 2 for d = 1.5s/, 
which is still acceptable for fabrication with electron beam lithography 2 and scanning probe 
lithography^ 2 -. 

VI. FDTD SIMULATION FOR NW WITH DIFFERENT SLOT WIDTHS w 

In our design, our aim is to set the optimum width of the air slot so that we can easily 
manipulate the NW in the center of the air slot without any loss of Q. Fig. [5h and b 
show the simulated electric field E x for different w values. The mode distibutions are very 
sensitive to w, especially from w = 50 nm to w = 75 nm. As w becomes larger, we observe 
that the electric field is enhanced in the air slot while the mode distributions become narrow 
and spread in the waveguide direction. This is similar to the case with different L values 
and an entirely drilled air slot described in Ref.—. The transversal profile of E x along the 
x-direction is shown in Fig. [5b. As with the entirely drilled air slot in Ref.—, the electric-field 
maxima are positioned at the boundary of the NW because the enhanced field originates 
from polarization charges induced in the NW walls. This effect is more pronounced when 
the slot walls have a boundary with the air of 50 nm< w ^ 150 nm. The intensities at the 
slot walls and in the center increase along with the increase in w. Here we simulated the 
square cross-section NW, but the trend may not be so different for the other type of NW 
facet. Fig. [5H shows a significant increase in Q from the design where 50 nm< w ^5 150 nm. 
The largest Q of 1.8 xlO 6 for w = 150 nm constitutes approximately a 6-fold improvement 
compared with that of the NW with a ID photonic crystal cavity proposed in Ref.—. For 
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w > 150 nm, Q decreases as a result of the decrease in the refractive index, i.e. the condition 
for total internal reflection (TIR) becomes severe as the slot width increases^. Additionally, 
we observe a blue shift of A and an increase in V c towards a large w value. The Q/V c derived 
from Fig. [SU is relatively constant for the design where w ranges from 75 nm to 150 nm. 
This means that we can choose the slot width within this range. However, another factor 
for decreasing Q, which limits the selection of w, is discussed in the following section. 

VII. FDTD SIMULATION FOR THE MISALIGNMENTS OF THE NW 
FROM THE CENTER OF THE SLOT 

We have already shown that a large Q can be obtained if the NW is combined with a 
large air slot. A large air slot has advantages in that it is easier to put NW into it. However, 
Q depends on the position of the NW (^-direction). This means that a large air slot, despite 
its high Q, may experience a large decrease in Q if the NW is not in the center of the air slot. 
Fig. [6^ and b show the calculated electric field E x for the NW inside 150- and 100-nm slots 
in the x — z and x — y planes, respectively. The spatial misalignment is indicated by p x . We 
observe that the mode distributions depend strongly on the position of the NW and they 
become broader when the NW is closer to the slot wall. A narrow slot width (w = 100 nm) 
provides stronger light confinement than a wider slot (w = 150 nm) for the same NW spatial 
misalignment (p x = —25 nm). The electric- field profiles for different spatial misalignments 
of w = 150 nm and w = 100 nm are shown in Fig. [6b and d, respectively. The maximum 
of the field shifts to the opposite direction of the NW misalignment and the field strengths 
exhibit monotonic fade out for larger NW misalignments from the center of the slot. The 
cavity wavelength (A) undergoes a red shift showing the potential for detuning using the 
NW position. When we compare the Q of the maximum misalignment (the NW attached 
to one side of the slot walls \p x \ = (w — si)/ 2) with that of the NW in the slot center, we 
find that the Q exhibits an approximately 7-fold decrease for w = 150 nm while that of 
w = 100 nm exhibits a loss factor of 2.4. From Fig. [SU, we found that the Q/V c values for 
the designs with w between 75 and 150 nm are almost the same, while the results for the 
misalignment of the NW in w = 100 nm and w = 150 nm shows that a narrow air slot is 
better for maintaining a high-Q in our design. Although we do not show the results for the 
misalignment of the NW for w values between 75 and 100 nm, we may conclude that the 
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optimum air-slot width for Q/V c can be found in this range. 

VIII. &-SPACE FIELD DISTRIBUTION OF THE CAVITY MODE FIELD 

Here we investigate the origin of high-Q with a large air slot and the decrease in Q with 
NW position misalignment through the x— and z— directional spatial Fourier transformation 
(FT) of the cavity field. Fig. [7J shows the FT spectra of the field profile in the middle of a 
photonic crystal slab for (a) si = w = 50 nm, (b) si = 50 nm, w = 150 nm, p x = nm, and 
(c) si = 50 nm, w = 150 nm, p x = —50 nm. The light cones of air are represented by white 
circles with radii of 2tt/X. Inside these circles, the TIR condition is broken for decomposed 
plane momentum components (k x ,k z ). As the slot width w increases, V c increases and the 
number of momentum components inside the light cone decreases (see Fig. [7Ji, b, and Fig. 
[5H). This condition results in the 4.5-fold improvement in Q. On the other hand, the NW 
misalignment causes an increase in the number of momentum components inside the light 
cone (see Fig. [7b). In this case, a 7-fold deterioration of Q is observed (see Fig. [5U and Fig. 

Et). 

IX. LIGHT-MATTER INTERACTION OF THE EMBEDDED EMITTER 
NW 

Here we focus on the light-matter interaction of the embedded emitter NW (QD NW) 
inside the slot of the line defect of the photonic crystals. There have been many studies of 
heterostructure QD NW, mostly concerning InAs/GaAs and InAs/InAsP QD NWs— ~— . 
Of these studies, two reports show that a photonic NW can control the spontaneous 
emission^ 4 . They found that the spontaneous emission was greatly inhibited (between 
12- and 16-fold) for a small-facet NW (si < 100 nm), from which they determined that 
the quantum efficiency of the QD NW was between 92 and 97 %22^2i. In terms of the very 
efficient funnelling of the spontaneous emission into the guided mode, this photonic NW 
is generally to be an efficient single photon source. However, the inhibition of the spon- 
taneous emission can be a drawback when using this small-facet photonic NW as a single 
photon source. Therefore, by using a photonic crystal cavity, we can expect the spontaneous 
emission rate of the small-facet NW to be enhanced. 



9 



For the following discussion, we will use the parameters of an InAs QD in InP NW since 
InAs/InP is already known to have potential as a single photon source in the telecom band^. 
To realize optimal coupling, we define the position of the QD as being the center of the NW 
while the transition dipole moment should be aligned with the electric field dipole. For a 
NW located in the center of the air slot, the center of the NW is also the location of the 
electric field maximum. 

As first noted by Purcell^, the radiative emission rate of an emitter coupled with the 
resonant cavity mode can be derived directly from Fermi's golden rule: 

2,71 f°° i \ 
T = ¥j (\Pa- ^e)\ 2 ) Pc(oo) Pe (uj)du (1) 

where p c {w) is the density of photon modes in the cavity, p e (w) is the mode density for 
the dipole transition, p a is the atomic dipole moment, and E(r e ) is the electric field at the 
location of the emitter normalized by a factor £ 2 = % r ^l,^ .» so that it gives to the zero 

2 J e^r)E I {r)d i r 

point energy. 

The spontaneous emission rate can be enhanced in two ways. The first is to increase 
the cavity mode density p c {oS), which is better known as Q. The other involves increasing 
the normalized electric field at the emitter position (£E(f e )), which is proportional to the 
decrease in V c . Therefore, the common figure of merit for the resonant cavities is the ratio 
Q/V c , which can be seen from the Purcell factor F p . In our case, we regard F p as the ratio 
of spontaneous emission rate of QD NW in the cavity ^"qd-nw compared with that of the 
QD in the bulk Tqe>; 



r r ,,„-X\l s , TQD-Nw(sl/Xem) ^ QD-NW , .-, , 

t v = b v x r = — p {*) 

1 QD 1 QD 

3QA 3 e(f )|i(f )| 2 



4vr% w 3 j e (r)\E{r)\ 2 d 3 r 
3QA 3 6Q 



(3) 
(4) 



4ir 2 n NW 3 V c tt 2 V; 

where fo is the center of the photonic crystal slab, which is the peak field and A is the cavity 
wavelength. pcav-NW j g ratio of the spontaneous emission rate of a QD NW in the 
cavity compared with that for a QD NW without a cavity Tqd^nw, which depends on the 
facet size normalized to the emission wavelength of the NW (s//A em ). In a small-facet NW 
(sl/X em < 0.16), the emission rate is inhibited (Xqd-nw I^qd < 1) while in a large-facet 
NW (sl/X em > 0.16), the emission rate is slightly enhanced (Tqd-nw I^qd ~ l}^ 1 ^. We 
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also define the dimensionless mode volume V c as 



( 




3 



v c = v c 



(5) 



Here we analyse the enhancement in the spontaneous emission rate caused by the increase 
in Q for the large slot width as previously seen in Fig. |5j Fig. [8] shows the dimensionless 
mode volume V c and the Purcell factor F p as functions of the position of the NW in the 
x— direction p x and the slot depth d for slots with w = 150 and 100 nm. For the slot with 
w = 150 nm, V c is slightly more sensitive than that in u> = 150 nm (see Fig. [8k and c). The 
minimum and the maximum V c values in the w = 150 nm slot are 5 and 10, respectively. 
Using a Q of between 3xl0 5 and 2.2xl0 6 and V c , we obtain an F p of between 1.5xl0 4 and 
2.1xl0 5 (see Fig. [8b and d for w = 150 and 100 nm, respectively). The position of the NW 
in the x— direction has a stronger influence on F p than the slot depth. For both large slots, 
the F p maxima are found for the NW in the center of the slot p x = nm, and the shallow 
depth slot d = 25 nm. The slot with w = 100 nm exhibits a more stable F p with p x than 
the slot with w = 150 nm. Here, we found the best F p of 2.1xl0 5 for w = 100 nm, p x = 
nm, and d = 25 nm, but we still may obtain a larger F p for 75 nm < w < 100 nm. 

As we discussed in relation to Eq. El F p in our calculation includes Fp av and Tqd-nw/^qd- 
This means that we can compare the emission rate enhancement for a QD NW before and 
after it is placed in the air-slot. Taking account si = 50 nm, A em = 1400-1600 nm, 
and assuming high quantum efficiency of the QD NW (> 90%), we can estimate that 
Fqd-nw/Fqd < 0.0625 22 . This estimation yields Fp av > 3.4xl0 6 and such an enhance- 
ment will greatly improve the single photon rate of the small size NW. Thanks to the 
large F p and Fp av , the present design provides a very good platform for cavity QED with 
a QD. We can expect a greater enhancement of the spontaneous emission rate for a QD at 
telecommunication wavelengths than that reported in Ref.— . 

The NW-photonic cavity may even enter the strong-coupling regime of light-matter in- 
teraction, where there is a coherent exchange of energy between the photon trapped in the 
cavity and the exciton trapped in the QD2&2S This occurs when the cavity field k, the QD 
exciton decay rate 7, the phonon dephasing rate 7^, and the cavity-QD frequency detuning 
S are smaller than the exciton-photon coupling strength g; 



K + 7 



+ ld P + i5 <2g 



(6) 



2 
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where g is simply written as 



qd / 3cA 2 



n NW V {2kT qd )AV c y ! 



Assuming that the radiative lifetime of the exciton of QD is 2 ns^ 4 - 1 ^ 6 -, which is converted 
into a radiative decay rate Tqd of 0.5 GHz, we obtain a g value of 170 GHz. The non- 
radiative recombination rate of the QD is about 42 MHz^2£, which is much lower than g. 
The dephasing rate 7^ for the telecom-band QD is typically 290 MHz^, and thus the only 
limiting factor shown in Eq. [6] is k. In our design, the smallest Q for the NW in the large 
slot is about 3xl0 5 . This is equal to k — 648 MHz, which is still much smaller than 170 
GHz. Therefore, we note that our design is suitable for the strong coupling regime even 
when the NW is misplaced in the slot of the photonic crystal. 



X. CONCLUSION 

We describe a novel photonic crystal cavity design, which is suitable for a position- 
controlled cavity for on-demand photonic controlled devices in photonic integrated circuits. 
This is the first design to realize a position-controlled cavity with a high Q of ~ 2xl0 6 and 
V c < 10. This design also shows that Q and V c are not sensitive to the fabrication of the slot 
(depth and width), the NW size, and an misalignment of the NW position. Our design is 
also useful for a cavity QED platform with an embedded QD NW. This system introduces 
the concept of the position-controlled cavity and emitter, which can be useful for coupling 
the emitter and the cavity. The NW position dependence of Q and the cavity wavelength 
are interesting as regards controlling the interaction between the cavity and the emitter. 
In terms of application, our system could be useful for the position-controlled devices such 
as position-controlled single photon sources, nanolasers, light emitting devices, and coupled 
cavities, though we need to do further investigation on the Q dependence with the spatial 
translation of the NW along the line defect. As regards fabrication, our approach requires 
only simple procedures involving the nanomanipulation of an NW inside a slot in a line 
defect of a photonic crystal. 
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Figure 1. (a) Schematic of a position-controlled ultrahigh Q-cavity using a single NW inside an air 
slot in a line defect of 2D photonic crystal. The inset shows the geometry of the sample. The white 
and the black letters indicate the parameters of the NW and the photonic crystal, respectively; a 
line defect width (W) of 0.98-2^/3, a hole radius (R) of 130 nm, a NW cross section side length (si) 
of 50 nm, an air slot depth (d) between 25 and 75 nm, a slot width (w) between 50 and 200 nm, 
and a NW length (L) between 3 and 5/^m. (b) Photonic bandstructure calculation for photonic 
crystals with a lattice constant (a) of 380 nm with only a line defect (black circles) and that with 
a NW of L = 3/itm and si = 50 nm and an air-slot of w = 50 nm and d = 50 nm (red squares). 
Other parameters are the same as in (a). A cavity formation created by introducing the air slot is 
shown by the green dashed line. 
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Figure 2. (a) Calculated electric-field (\E\) and magnetic-field (|-ff|) distributions in the x — z plane 
for a = 380 nm, L = 3 /urn, R = 130 nm, and si = w = d = 50 nm. The color scale is linear and 
blue filled and empty rectangles in the middle represent the NW and the air slot, respectively. 
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Figure 3. (a) Calculated z-component of the magnetic- field (|-f^ 2 |) distributions in the y — z plane 
for different L values of the NW and other parameters; a = 380 nm, R = 130 nm, and si = w = 
d = 50 nm. The color scale is logarithmic and photonic crystal slabs are depicted by solid lines. 
The tick lines show the lengths of the NW and the relative positions of the NW in the slab. Each 
field is normalized by the total field energy in the whole calculation area, (b) The magnetic-field 
profile along the z direction (the maximum field) indicated by dotted lines in (a), (c) Summarized 
performance of Q and V c for different L values. 
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Figure 4. (a) Calculated y-component of the magnetic-field fl-Hyl) distributions in the x — y plane 
for different d values of the air slot and other parameters; a = 380 nm, R = 130 nm, L = 3 //m, 
and si = w = 50 nm. The color scale is logarithmic and photonic crystal slabs are depicted by 
solid lines. Each field is normalized by the total field energy in the whole calculation area. The 
magnetic-field profiles along (b) the y-direction and (c) the x-direction correspond to the vertical 
dashed and horizontal dotted lines, respectively, indicated in (a), (d) d dependences of cavity peak 
A, Q, and V c . 
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Figure 5. (a) Calculated x-component of the electric-field (l-E^I) distributions in (a) the x-z and 
(b) x-y planes for different w values of the air slot and other parameters; a = 380 nm, R = 130 nm, 
L = 3 /im, and si = d = 50 nm. The color scale is logarithmic and photonic crystal structures are 
depicted by solid lines. Each field is normalized by the total field energy in the whole calculation 
area, (c) The electric-field profile along the x-direction corresponds to the vertical dotted lines in 
(a), (d) Cavity peak A, Q, and V c as a function of w. 
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Figure 6. (a) Calculated ^-component of the electric field (|-Er|) distributions in (a) the x-z and 
(b) x-y planes for different w values with a misalignment of the NW p x and other parameters; a 
= 380 nm, R = 130 nm, L = 3 /im, and si = d = 50 nm. The color scale is logarithmic and 
photonic crystal structures are depicted by solid lines. Each field is normalized by the total field 
energy in the whole calculation area. The electric-field profile along the x-direction corresponds to 
the horizontal dotted lines shown in (a) for (c) w = 150 nm and (d) w = 100 nm. Data for p x = 
were obtained from Fig. [St and the grey lines show the air slot dimension. 
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Figure 7. Spatial Fourier transformation spectra obtained from the x-component of the electric 
field (E x ) in the x-z plane for (a) a NW in a width-matched slot (si = w = 50 nm), (b) a NW 
in the center of a wide slot (w = 150 nm), and (c) a NW attached to one side of the slot walls 
(w = 150 nm, p x = —50 nm). The leaky regions are shown by solid white lines. Other parameters 
used in the calculation are a = 380 nm, R = 130 nm, L = 3 ^m, and d = 50 nm. 
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Figure 8. (a) Dimensionless mode volume V c and (b) Purcell factor F p for a slot width w of 150 nm 
calculated for different NW positions p x and slot depths d. (c) and (d) are V c and F p for w = 100 
nm, respectively. Other parameters used in the calculation are a = 380 nm, R = 130 nm, L = 3 
^m, and si = 50 nm. 
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